Astronomy & Astrophysics manuscript no. 2Shears printer 


© ESO 2008 


February 5, 2008 





O 
O 

(N 

c3 



Self-consistent theory of turbulent transport in the solar 

tachocline 

II. Tachocline confinement 

Nicolas Leprovost and Eun-jin Kim 
Department of Applied Mathematics, University of Sheffield, Sheffield S3 7RH, UK 



> 

(N 

in 

o : 
\o , 
o ■ 

^: 

Oh; 
6 1 



Received / Accepted 



ABSTRACT 



Aims. We provide a consistent theory of the tachocline confinement (or anisotropic momentum transport) within an hydrodynamical turbulence 
model. The goal is to explain helioseismological data, which show that the solar tachocline thickness is at most 5% of the solar radius, despite 
the fact that, due to radiative spreading, this transition layer should have thickened to a much more significant value during the sun's evolution. 
Methods. Starting from the first principle with the physically plausible assumption that turbulence is driven externally (e.g. by plumes 
penetrating from the convection zone), we derive turbulent (eddy) viscosity in the radial (vertical) and azimuthal (horizontal) directions by 
incorporating the crucial effects of shearing due to radial and latitudinal differential rotations in the tachocline. 

Results. We show that the simultaneous presence of both shears induces effectively a much more efficient momentum transport in the 
horizontal plane than in the radial direction. In particular, in the case of strong radial turbulence (driven by overshooting plumes from the 
convection zone), the ratio of the radial to horizontal eddy viscosity is proportional to JT i/3 , where is the strength of the shear due to radial 
differential rotation. In comparison, in the case of horizontally driven turbulence, this ratio becomes of order -e 2 , with negative radial eddy 
viscosity. Here, e (<k 1) is the ratio of the radial to latitudinal shear. The resulting anisotropy in momentum transport could thus be sufficiently 
strong to operate as a mechanism for the tachocline confinement against spreading. 
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1. Introduction 

One of the outstanding problems in solar physics is to un- 
derstand the dynamics of the tachocline, a thin layer where 
the transition from latitudinal rotation in the convective zone 
to uniform rotation in the radiative interior takes place. This 
is particularly true since the tachocline links two regions of 
very different transport properties, thereby playing a crucial 
role in the overall angular momentum transport and chemical 
mixing on the course of the solar evolution. The tachocline 
is also crucial for the solar dynamo since the strong shear 
in the tachocline is believed to take part in the process by 
whic h the magnetic field of th e sun is created. Helioseismic 
data JCharbonneau et ai1ll999h constrain the thickness of this 
layer to be only a few percent of the solar radius. This poses a 
challenging problem since a radially localized shear flow (as- 
sociated with the radial differential rotation) naturally tends 
to spread during the c ourse of solar evolution, as shown by 
LSpiegel & Zahnl i 19921) . and should have reached a thickness 
much broader than what is observed today. 

Some physical processes have been pro posed to be re- 
sponsible for the tachocline confinement. First, Spieg eT& Zahnl 



(1992) have shown that if the turbulence is highly anisotropic 
(being much more vigorous in the horizontal direction than in 
the radial one), the spreading of the solar tachocline could be 
limited. More precisely, it is the turbulent transport of angu- 
lar momentum which has to be anisotropic, meaning that the 
horizontal turbulent viscosity has to be much larger than the 
radial one. Physically, this would enable the differential ro- 
tation to be smoothed out by transport in the horizontal di- 
rection rather than in the vertical one. An alternative mecha- 
nism for the st abilizati on of the tachocline is due to magnetic 
field jRiidiger & Kichatin"ovlll997t iGough & McIntvTdTl998l: 
iMacGregor&ChaAomiemJ 1 9991) : even a rather weak poloidal 
magnetic field (of the order 10~ 4 G) in the solar interior could 
prevent the downward flow from the convection zone (driven, 
for example, by rotation-induced meridional circulation) from 
penetrating deeply into the radiative interior, thereby confining 
the differential rotation to a thin region of space and rendering 
a uniformly rotating interior. 

The relevance of these models (and any other that can be 
constructed with a physical process) to the sun is however not 
well asserted. In particular, the success of these models de- 
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pends crucially on the values of transport coefficients such as 
turbulent viscosities and magnetic diffusivities, which are of- 
ten crudely parameterized, unless their molecular values are 
used. Furthermore, the Gough & Mclntyre model is far to com- 
plex to be analyzed thorough ly and on ly a simplified version of 
it is actually tractable Jc.araudl EoO.l).' In the Spiegel & Zahn 
scenario, it is of prime importance to identify physical mecha- 
nisms which can lead to anisotropic momentum transport and 
then to derive eddy viscosities from first principles (or from 
Navier-Stokes equation). Note that while Spiegel & Zahn in- 
voked a strong stable stratification as a source of anisotropic 
turbulence, they did not derive the values of eddy viscosities. 

The purpose of this paper is to provide a consistent theory 
of the anisotropic momentum transport. The physical mecha- 
nism that we invoke is the turbulence r egulation by shearing 
effect , the so-cal led shear stabilization JBurrelllfl997t |tlahm 
2002; iKiml 12004). It is basically because a shear flow acting 
on a turbulent eddy creates small scales in the direction or- 
thogonal to the flow and thus enhanc es dissipation, reducing 
the transport and turbulence intensity (Kim&Dubrulle 2001; 
Kim & Diamond!l2003t iKim et aljEoo4 lKimll2006l) . Indeed. 
KimI 2005) has shown that a shear due to a stable radial dif- 
ferential rotation in the tachocline not only suppresses turbu- 
lent transport and turbulent intensity, but also leads to an "ef- 
fectively" anisotropic transport of particles via stronger reduc- 
tion in radial transport compared to horizontal one. However, 
the key question of anisotropic momentum transport, n eces- 
sary fo r the tachocline confinement, was not addressed in lKirrJ 
ll2005l) since the effect of latitudinal differential rotation in the 
tachocline was simply neglected for simplicity. In this paper, 
we provide a theory of anisotropic momentum transport in the 
tachocline within an hydrodynamical turbulence model, by tak- 
ing into account the crucial effects of both radial and latitudinal 
differential rotations in the tachocline. Specifically, we com- 
pute the transport properties of turbulence self-consistently un- 
der the physically plausible assumption that turbulence arises 
by an external forcing (e.g. from plumes penetrating from the 
convection zone). Most results will be derived in the limit of 
strong shear, but with latitudinal she ar wea ker than radial shear, 
as relevant to the tachocline. As in lKinj J2005I) . we shall ne- 
glect the effects of stratification, magnetic field and rotation to 
elucidate the crucial role of shearing due to radial and latitudi- 
nal differential rotations in momentum transport. In particular, 
we shall show that shear alone can lead to anisotropic momen- 
tum transport, with a more efficient horizontal transport. While 
stratification, magnetic field and rotation can also contribute to 
anisotropic turbulence, these are outside the scope of this study 
and will be addressed in future papers. 

The remainder of the paper is organized as follows: we pro- 
vide our model of the tachocline and solve the governing equa- 
tions in section|2j then study the turbulence amplitude and eddy 
viscosities in section[3]and section |4] respectively. Sectionals 
devoted to our conclusions and discussions on the implications 
for the solar tachocline dynamics. Appendices contain some 
details of the algebra involved in solving our model and calcu- 
lating the turbulent amplitude and viscosity. 



2. Model 

To elucidate the effect of radial and latitudinal differential ro- 
tation on turbulent transport, w e adopt a simplified model for 
the tachocline as in Em|(2005) and use a local Cartesian refer- 
ence frame where x, y and z denote local radial, azimuthal, and 
latitudinal directions, respectively. We capture the radial and 
latitudinal differential rotation by a large-scale velocity field 
given by V (x,z) = -{xJK x + zM z )y. Here, M x = d x U > 
and M z = z Uq > are the strength of the radial and lati- 
tudinal shear due to radial and latitudinal differenti al rotation , 
respectively. In the quasi-linear approximation (Moffatt 1978), 
the Navier-Stokes equation for the small-scale fluctuating ve- 
locity can be written: 

d,\ - (x$l x + zJi z )d y \ - (M x v x + W z v z )y = -Vp + vV 2 v + f , 

V- v = . (1) 

Here f is the forcing at small scale and v is the viscosity of 
the fluid. We will use a Fourier-transform with time-dependent 
wave number to accou nt for the ef fect of shearing on eddies 
non-perturbatively re.g. lKiml J2005l) l: 

Here, both radial and latitudinal wave-numbers evolve in time 
as follows: 

k x (t) = k x (0) + M x t and k z (t) = k z (0) + M z t . (3) 

To absorb the viscosity term, we set V = Vexp[v(2(f)] where 
2(0 = [k^/lkyMx+k^t+kl/Sky&lJ. The quasi-linear equations 
for the fluctuating part of the velocity field then become: 

d,v x = -ik x p + f x , (4) 

8,Vy - (Ji X V x + Ji Z V Z ) = -ikyP + fy , 

d,v z = -ik z p + f z , 

= k x v x + kyVy + k z v z . 

To solve equation @, we change the time variable from t to 
t = k x {t)/k y and introduce new variables: J[ = ${ x , e = M z /A x 
and <p = [k z (Q) - ek x {Q)]/k y . For parameter values typical of the 
solar interior, we have: 31 = AUq/AX ~ A£l/(h/Ro) = 3 x 
10~ 6 s _1 for the tachocline of thickness 4% of the solar radius, 
and M z = AU/Az ~ 4AO ~ 8 x lO^s" 1 . With these values, 
the ratio of the azimuthal and radial shear is a small parameter 
e ~ 2.7 x 10~ 2 . We shall thus calculate the turbulent amplitude 
(section and eddy viscosities (section |4} up to order e 2 . In 
terms of the new variables, Eq. (0} becomes: 

Wd T v x = -ikyTp + f x , (5) 

J[d T v y - J[(v x + ev z ) = -ikyp + f y , 

Jld T v z = -ikyier + (f>)p + f z , 

= tv x + Vy + (er + <p)v z . 

The solution of this system of equations can be found after a 
long but straightforward algebra with the result (see appendix 
EJfor details): 

r , £i(Ti) f y + e<pT e<p r -n 

J To M {(y + e 2 )R(T) (y + e 2 ) 3 / 2 L J J 
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- I dr\ 
= I dr\ 

Jt„ 



h 2 (ri) 



31 y+e 2 



(6) (v 2 z ) = J^^f d 3 k@ {^((g + c?) 2 4 + ga 2 }i!,u(k) 



(y + e 2 )R(r) I J T0 



Ji y+e 



{* d fii(Ti) < e-<f>T « 

J T0 Tl M \{y + e 2 )R{T) (y + , 

+ f rfr, 

Jt 



^1 l(y + e 2 )/?(r) (y + e 2 ) 3 / 2 

4(T!) 1 



t(t) - r(n) 



y + e 2 ' 
Here, 

fl( T ) = (e 2 + l)r 2 + 2e^r + y , y = 1 + 2 , (7) 

ty n t ( (e 2 + l)T+e<p \ 

T(t) = arctan — , 

h(T) = [y + 0er]/, - [(1 + e 2 )r + 0e]£ + (e - 0t)£ , 

4(t) = -e/»(T)-0/,(T) + / z (T). 



3. Turbulence amplitude 

We first examine how turbulence amplit ude is affecte d by lati- 
tudinal shear, comparing the results with lKiml J2Q05). In order 
to compute (v 2 ), we assume the forcing to be incompressible 
(thus all the quantities can be expressed in terms of the x and 
Z components only) with the statistics that is spatially homoge- 
neous and temporally short correlated with the correlation time 
Tf. Specifically, we take: 

</Kki, tdfjOti, f 2 )> = Tf (2n) 3 6(k 1 + ka) m - h) ^y(ka) , (8) 

for i and j = 1 , 2 or 3. The ifr^ functions are the power spectrum 
of the forcing. 

A typical solar value of the molecular viscosity is v ~ 
10 2 cm 2 s _1 and thus the parameter ^ = vk 2 /Jl is a small quan- 
tity provided that k y < 10~ 4 cm ~ 10~ 6 /#o where Ho is the 
pressure scale height at the bottom of the convection zone. The 
relevant length scale in the azimuthal direction being proba- 
bly larger than 10 4 cm, this permits us to consider the strong 
shear limit £ <K 1 . In addition, since e is a small parameter, as 
noted previously, we expand the result up to e 2 , keeping only 
the dominant term in £ and e for each component of the forc- 
ing ftij. A long but straightforward algebra, very similar to that 
of the calculation of the turbulent viscosities (see next section 
and appendix[E}, then gives us the turbulence amplitude in the 
strong shear limit £ <s 1 to leading order in ^ and e: 



2ba 

+ ^1300 + #33(k) 



Here, 



k x /k y , b = kjky , g = 1 + b 
n/2 - arctan(a/ yfg) = arctan( yfg/a) , 



(10) 



0o = r(l/3)/3(3/2£) 



1/3 



f = Vkt/JTI . 



In Eq. (|9), we kept only the leading order contributions in ^ and 
e for the terms proportional to i^n, 1^13 and 1^33. For instance, 
for the component proportional to 1^33 in the radial turbulence 
amplitude, we expanded the quantities up to order 0(e 2 ). It is 
worth noting that (v 2 ) is always positive, as it should: while the 
second term proportional to tfr^ is not always positive, in order 
for the sum of the two terms proportional to and 1^13 to be 
negative, t/03 has to be at least e _1 larger than (Ai 1 ; in that case, 
0-33 would be at least of order e~ 2 and thus the last term would 
dominate the other two terms, making (v 2 ) positive (recall that 
£ <sc 1 and thus @q » 1). In comparison, the horizontal turbu- 
lent velocities in Eq. © are obviously positive and are the same 
as that obtained by Kim (2005) in the case without a latitudinal 
shear. 

Eq. indicates that the amplitude of the three compo- 
nents of the velocity depends not only on the radial and lati- 
tudinal shear (shearing rate J[ and e$[) but also on the typical 
wavenumber k and the power spectrum of the forcing 0v/(k). 
We first examine how various components are affected by the 
radial and latitudinal shear. As can easily been seen from Eq. 
(|9), the amplitude of all three components of the turbulent ve- 
locity is reduced due to radial shear, becoming very small as J[ 
increases. Specifically, the horizontal turbulence is reduced by 
a factor J[Qq 1 °c ^t 2 ^ 3 , the inclusion of the latitudinal shear 
having no effect on the horizontal turbulence amplitude. In 
comparison, the radial turbulence velocity has a small contri- 
bution of order 0(e 2 ) from t/133 because of a latitudinal shear. In 
order to understand this, it is instructive to consider a new ref- 
erence frame such that the shear flow Uoy, depending on both x 
and z, becomes a function of only one coordinate xf . This coor- 
dinate transformation can easily be obtained by rotating x and 
Z axis by the angle 8 = arctan^j/y^) ~ e around the y axis 
(see figure^ as follows: 



cos 6 sin 6 \ I x 
- sin# cos 9 J \ z 



1 e 
-e 1 



(11) 



+ a n 



V? 



(9) 



(g + a 2 )e r a 1 -. 
H \-a + 



VI 



■/C 



J d 3 kg ^(ig + a 2 ) 2 4 + ga 2 y n (k) 



2b^a r, 
+ -Ai3(k) + fcV 33 (k) 



where we used e <k 1 . In the new (x',z') coordinates, the large- 
scale velocity field is given by = -ijX^+J^y^x'y', depend- 
ing only on the radial direction, thus becoming a purely radial 
shear. Thus, in the (x',z') coordinate, we can use the results of 
( MS (120051) that <vf > oc &T l g ~ yr 2/3 and <v' 2 > oc yr 1 . By 
] ( Ao(k) + e 2 Qo (A33(k) \ , relating these to (v 2 ) and (v 2 ) via Eq. (II It . we easily obtain the 
term proportional to 1/133 in (v 2 ) in Eq. (|9j- That is, the latter 
is the result of the dependence of (v' 2 ) in the (x',z') reference 
frame with an additional factor of e 2 ~ sin 2 8 coming from 
the projection of this term onto the x axis. We also note that 
(v' 2 ) - (v?) for f « 1, consistent with Eq. (|9}. We will further 
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arctan(jp) ~ e 



Fig. 1. Sketch of the coordinate transformation from (x, z) to 
(x',z')- In the new referential, the shear is purely radial (see 
main text for details). 

use this geometrical argument when discussing the turbulent 
viscosities. 

Keeping these in mind, we now discuss the influence of 
the different types of forcing on the turbulent intensity. If the 
turbulent motions are primarily driven in the radial direction 
(0-33 = (^13 = 0), the latitudinal shear has no effect, and the re- 
sults are the same as in the case of a radial shear only: the radial 
turbulence is reduced by a factor J[ and thus, the ratio of the ra- 
dial to horizontal turbulence amplitude (taking the y direction 
as well as the z direction) is of order '. Alternatively, if 
the small-scale flow is mainly driven in the horizontal direction 
{\fi\i = if/13 = 0), this component {v 2 x ) is reduced by a factor 
of e 2 yr 2/3 , with the ratio (v 2 )/(v 2 ) oc e 2 , for the reason that 
was just discussed. Again, note that this is different from the 
case without latitudinal shear where the radial turbulence in- 
tensity was null (at the order of the approximation). The ratio 
(v 2 )/(v 2 ) oc e 2 is still small in the case of the sun but interest- 
ingly depends only on the ratio of the strength of the azimuthal 
and radial shear but not on their absolute magnitudes. 

4. Turbulent viscosity 

We now proceed to the calculation of the radial v" and hori- 



zontal viscosity vj defined as: (v x v y ) 



x d x U = vfM and 
(v z v y ) = -v^d z Uo = vfe^l. From the previous section, one 
might naively estimate the ratio of these two to be of order: 



(v*Vy) 

(VyV z ) 



M> ,1,6 

\(vl) ~* 



(12) 



for a jc-dominant and z-dominant forcing respectively. 
However, the momentum fluxes (v x v y ) and (v y v z ) depend not 
only on the velocity amplitudes, but also how closely differ- 
ent components of the velocity are correlated (i.e. the phase- 
relation). That is, the ratio of the two should be given as fol- 
lows: 



(VxVy) _ | (V 2 ) COS 6 xy 

j(v 2 ) cos,6 zy 



(VyV z ) 



(13) 



where cos 6 xy and cos 6 zy (the so-called cross phase) represent 
the phase-relation between x and y and z and y components of 
the velocity, respectively. Thus, the estimate given in Eq. fl!2i 
may not be true if the cross-phase cos 6 xy and cos 6 zy are af- 
fected in an anisotropic way due to the different strength in the 
radial and latitudinal shear. We shall show that this is indeed the 
case. Specifically, cosS xy / cos<5 zy oc g 1 / 6 when the turbulence 
is driven mainly in the x-direction while cos<5 xv / cos6 zy = -1 
when driven in the horizontal direction. 

After a long, but straightforward algebra, we can obtain 
eddy viscosities in the following form (see appendix|B]for de- 
tails): 



4eb 3 a 

"~vf 



■KOI 



£o]<An(k) 



(g + a 2 )2eb 2 2e 2 b 2 (g + a 2 ) _ . 



,3/2 



8 

T f C , 1 I g + a 2 r 4eb 4 n 
(2^J rfk e{V^ + ^°^ Il(k) 



2b 2 (g + a 2 ) 



7 3/2 



^ 13 (k) + efe(k)) 



: k x /k y , b — k z /k y , 
arctan( yfg/a), g 



= 1 + b 2 , KQ = tt/2 - 
r(l/3)/3(3/2£) 1/3 and 



Here, again, a 
arctan(a/ ^[g) - 
% = Dtf/M. 

Eq. dl4> clearly shows that both v" and vf are reduced by 
shearing, becoming very small as J[ increases. The exact scal- 
ings of the eddy viscosities however depend on the properties 
of the forcing since Qq oc ^'^ 3 oc appears in some of 

the coefficients proportional to </<n and if/^ in Eq. (I14> with 
both positive and negative signs, making the estimate rather 
complex in the case where the forcing is mainly driven radially 
by iff ii. In order to obtain a transparent scaling relation in this 
case, it is illuminating to consider a physically plausible case 
where the turbulent plumes, coming from the convection zone, 
are highly elongated in the x-direction, with small-scale struc- 
tures in the radial direction being mainly created by radial shear 
(see figure|2)- In this case, it is very likely that the forcing will 
select small wave numbers in the radial direction, permitting 
us to crudely set a - k x /k y ~ in Eq. (I14> . Mathematically, 
this is justified if the support of the iffy functions is localized 
near the line k x /k y = 0. In that case, the radial viscosity does 
not contain terms involving e or £ while the horizontal compo- 
nent is proportional to go = r(l/3)/3(3/2£)^ 3 . Therefore, the 
ratio of these two viscosities is of order v^/vf oc ^^ 3 «; 1. 
In other words, the momentum transport becomes anisotropic 
with a much larger horizontal transport compared to the radial 
one. 

Interestingly, the contribution from 1^33 to v x T x comes with 
the opposite sign to that from \fi\\ while the contribution from 
i/<33 to come with the same sign as that from ip\\. Thus, as 
0-33 takes a non-vanishing value due to horizontal forcing, it 
will further decrease while increasing until v£ l vanishes 
and then becomes negative. In other words, the anisotropy in 
momentum transport becomes even stronger with v x T x jv^ <K 
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Fig. 2. Effect of the large-scale shear (left) on plumes com- 
ing from the convection zone (right). In addition to the tilting 
induced by the large-scale shear, the turbulent velocity field in- 
duces stochastic distortion of the eddies which is not drawn 
here for simplicity. 

f 1 ^ 3 . In the extreme limit where the turbulence is solely driven 
horizontally with — if/^ — 0, v X j jv-j oc -e 2 . Note that this 
is true in general, regardless of the form of power spectrum of 
(A33. Since e 2 <K 1, the ratio of the magnitude of the two is 
small. It is important to note that the two viscosities have dif- 
ferent signs, with the radial (vertical) viscosity being negative 
while the horizontal one is positive. To understand these re- 
sults, it is useful to consider the coordinate (x',z') given in Eq. 
dl 11 1 where the shear flow depends only on x' . By an elemen- 
tary calculation, we can easily show that the turbulent viscosity 
Vj X in the (x',z') coordinate can be expressed in terms of v" 
and vip in the (x, z) coordinate as follows: 

^ = ^fT^f[^ + ^ 2v? ] ■ (15) 

We now recall that in the (x',z') coordinate with a purely ra- 
dial shear, v x T x — in the c ase of a pure ly horizontal forcing 
(i//n = if/n - 0), as shown by Kim (2005). Thus, requiring that 
the right hand side of equation (I15> be zero immediately gives 
us the result v^/vf oc -e 2 , as found in this paper. This is an 
interesting result since the negative viscosity here is not related 
to the 2-dimensionality of the flow in any sense, being merely 
due to the fact that a special rotation of the Cartesian frame 
maps our problem to the case of a purely radial shear. 

To summarize, for a reasonable choice of the forcing, 
Vj jv-j oc ^ 3 <g 1 and -e 2 for radial and horizontal forc- 
ing, respectively. In these two limiting cases, the ratio of the 
cross-phase given in Eq. (1131 satisfies cos 5 xy l cos 5 zy oc £ x / 6 
and -1, respectively. It is worth noting that the small value 
of cosSxyl cos (5^ oc £'/ 6 in the former case physically makes 
sense since a stronger decorrelation between velocity compo- 
nents in the radial direction and horizontal plane is expected 
because of a stronger radial shear. These results thus suggest 
that horizontal momentum transport can be much more efficient 
than radial transport due to radial and latitudinal differential ro- 
tations. 



5. Conclusions and discussion 

We presented a novel mechanism for the anisotropic angu- 
lar momentum transport to explain tachocline confinement. 
Starting from the first principle, we have consistently derived 
eddy viscosities and turbulence amplitude by incorporating the 
crucial effects of both radial and latitudinal differential ro- 
tations. We have shown that even in the absence of stratifi- 
cation and magnetic field, the shear induced by differential 
rotation alone can lead to anisotropic momentum transport. 
Specifically, when the turbulence is mainly driven radially by 
overshooting plumes coming from the convection zone, the ra- 
tio of the radial to horizontal eddy viscosities can be propor- 
tional to ^ 3 , becoming very small as the radial shear J{ in- 
creases. Here, £ = vk 2 /^l and J[ is the shearing rate in the ra- 
dial direction. In comparison, in the case where the turbulence 
is mainly driven horizontally, this ratio is proportional to -e 2 
where e («: 1) is the ratio of the shearing rate in the latitudinal 
direction to that in the radial one. In this case, the radial eddy 
viscosity becomes negative. These results thus suggest that hor- 
izontal momentum transport can be much more efficient than 
radial transport, with a strong anisotropic momentum transport, 
as a result of the shearing effects due to radial and latitudinal 
differential rotations. 

As mentioned in introduction, the anisotropic momentum 
transport has important implications for the dynamics of the 
tachocline. Spiegel & Zahn ( 1992) have shown that the radia- 
tive spreading of the tachocline could be limited if the ratio of 
the radial to horizontal diffusivities satisfies the following rela- 
tion: 

V T lh\ 2 , 

4r <K — ~ 3x 10~ 3 , (16) 
vf \r J 

where ro is the radius at which the tachocline is located and 
h its thickness. The estimate on the right hand side has been 
obtain ed using helioseismologic values of ICharbonneau et alJ 
(1999). Indeed, this criterion can easily be satisfied for a rea- 
sonable choice of the forcing as follows. First, in the case of 
radial forcing which is dominated by the components with 
k x lky <K 1, the Spiegel & Zahn criterion is easily met for 
a characteristic length-scale in the azimuthal direction L v = 
Inky 1 » 2Ti(r /hf(v/Jiy 12 ~ 3 x 10 5 m for typical values 
of the solar tachocline J[ ~ 3 x 10~ 6 s _1 . In the case of hori- 
zontal forcing, the ratio of radial to horizontal eddy viscosities 
becomes of order e 2 ~ 7 x 10~ 4 . Furthermore, in this case, the 
radial eddy viscosity becomes negative, reinforcing the pres- 
ence of radial shear with sharp gradient in a localized region. 
Therefore, these results suggest that the anisotropy in momen- 
tum transport could be strong enough to operate as a mech- 
anism for the tachocline confinement against spreading. Our 
results, however, cannot offer a mechanism for a uniform 
rotation in the solar interior. As noted in the introduction, 
this could be due to the presence of magnetic fields in this 
region (Riidiser & Kichatinov 1997; Gough & Mclntvre 
1998; MacGregor & Charb onneaull999h . 

We note that in this paper we have neglected some effects 
which could be important in the lower part of the tachocline. 
For example, the presence of stratification or rotation can gen- 
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erate a wave-like turbulence and consequently severely affects 
its transport properties. The effect of Coriolis forces has been 
studied in the context of the transport of angular momentum in 
the convection zone to elucidate the permanence of a differen- 
tial rotation profile in that region. It has been shown that, in the 
strong rotation limit, the turbulent viscosity is more reduced by 
a factor of 4 ( Kichatin ov et all 1 994) in the direction perpendic- 
ular to the rotation axis compared to the one parallel to the rota- 
tion axis. Furthermore, rotation induces a non-diffusive part in 
the turbulent viscosity (the so-called A effect) which prevents 
the solid rotation to be a solution of the turbulent momentum 
equation. We can expect a density stratification to have a sim- 
ilar ef fect to the Coriol is force as the linear equations are the 
same (Greenspan 1968) once the vertical and horizontal direc- 
tions are interchanged. In that case, we expect the component in 
the direction of the stratification to be more reduced than in the 
plane of constant density. In fact, such an anisotropic eddy 
viscosity has been observed in a recent numerical simula- 
tion of a stably stratified turbulence driven by pe netrative 
convection with an imposed shear (Miesch 2003). In par- 
ticular, the angular momentum transport was found to be 
diffusive in the latitudinal direction while anti-diffusive in 
the radial direction, with turbulence mixing potential vor- 
ticity [see also Mclntyre (2003)]. Another important effect is 
due to magnetic fields in the tachocline (as suggested by solar- 
dynamo models) by imposing Alfven-like structure in the prop- 
erties of turbulence. Furthermore, the simultaneous presence 
of toroidal magnetic field and latitudinal differential rotation 
can lead to a large- scale joint instability ( Oilm an & Foxll 99*71 
iDikpati et all2004") which could further enforce an anisotropic 
turbulence. The inter play between shear and r otation is under 
current investigation (Lepr ovost & KinjEo0 6) and the related 
problems of stratification and magnetic field will be published 
in future papers. 

Appendix A: Solution of system © 

Using the new field u — v x + ev z , we obtain the following equa- 
tions: 

3id T [R(T)u] = In , (A.l) 
tf[y + 4>eT]d T v x = -nd T [T 2 u\ + [y + 6er]f x - r/ v - </rrf t , 

where y = I + (f> 2 , R(f) = (e 2 + 1)t 2 + letpT + y and ti\ - 
[y + <per]f x - [(1 + e 2 )r + <pe]f y + (e - <frf) f z . The first equation 
can be readily integrated to obtain u — J dr\ Ii\(t\)I3UI(t\). 
Then, the second equation can be used to obtain v x . Here, we 
provide some of the main steps in the calculation. First, we 
calculate the first term on the RHS of Eq. (IA.1> : 



(A.2) 



2 . , _ r 2 ln (t) 2r[y + ecf>T] f T h\ (n ) 



d T [r l u] = 



t 2 /h(t) 2r[y + 6^T] f T 

am + m 2 J Tn 



Then plugging Eq. iA.H in the second equation of JA. li . we 
obtain: 



Jr„ R(T2) 2 Jt 

+ f 



= _ r dTl Mm r 

•J To 



dr 2 



2t, 



R(r 2 ) 2 



[R(n) - j\\f x - nf y - Ti(# + en)/, 



&R(t x ) 



J T0 Tl M \(y + e 2 )R(T 

•JTn 



m L(y 
hi{T\) e 
a y + e 2 



)R(t) (y + e 2 ) 3 / 2 



T(t) - T( Tl ) 



Here we used the formu la (2.175) and (2.172) of 
iGradshtevn & Rvzhikl i 1965b to calculate the inte- 
grals; T(j) - arctan([(e 2 + l)r + etp]l ^y + e 2 ) and 
/i2(t) = —ef x (j) - <f)f y (j) + j z (t). The two other com- 
ponents of the velocity can easily be obtained by using 
v z = (u - v x )/e and v y = -tv x - (er + <p)v z . 

Appendix B: Calculation of the turbulent 
viscosities 

We symbolically write (v x v y ) and (v z v y ) as follows: 

< V ^v> = TTTT^ f d 3 k{Vll<l>U + Vl20l2 + V22022) , (B.l) 



(v 



<- v v> = ,n, T L m \ d 3 k{H u 4> u + H n <pi2 + H22022} 
(2n) s M J 



where the 0y functions are the power spectra associated with 
the forcing functions hi and h^. 

Mk^hjik^ti)) = r / (27r) 3 5(ki+k 2 )5(f 1 -f 2 )0, 7 (k 2 ) ,(B.2) 

for i and j = 1,2. The functions Vy and ffy (for i and j = 1,2) 
are expressed as integrals over the variable t, for example, Vn 
can be written: 



Vv 



•Ju 



^ Te -2v{e(T)-fi(a)}| 

2e<p 2 



(e 2 - y)4> + e[e 2 + 1 



2 \r 



(y + e 2 ) 



2\2 



(y + e 2 ) 5 ' 2 



T(t) - T(a) 



(B.3) 



where, a = fc r /fc v , <p = (k z - ek x )/k v and y = 1 + tp 2 . Following 
iKinj (12005b . this integral will be estimated in the strong shear 
limit, i.e. for f = vk 2 /Jl <K 1. Thus, v{Q(t) - Q(a)} = £[(e 2 + 

l)T 3 /3 + e0r 2 +yT-{ya + (e 2 +l)fl 3 /3+^ra 2 )}] where b = k z /k y . 
Taking the limit £ — > 0, we evaluate each term in (IB.3I by 
keeping only terms to leading order in The results are: 



3\ 



(A3) 



[R( Tl ) - r\}f x - Ti/ y - ti(0 + en)f y 

^«(Tl) 



Via 
V 22 



y + £(f>a 



e<t> 



2(y + e 2 ) 2 R(a) 2(y + e 2 ) 5 / 2 2(y + e 2 ) 



(y + e 2 ) 3 
(e 2 - y)0 
(y + e 2 ) 5 ' 2 



K 2 + 



-l-e 2 ) 



(y + e 2 ) 5 /2( e 2 + 1) 
2e<^ 2 e(l+e 2 -4> 2 ) 



• K0 + 



( r + e 2)5/2 - ( r + e2)2 (e 2 + 1) 



(y + e 2 ) 2 
6a — e 



(B.4) 







K + 



eS 2 



2(y + e 2 ) 2 R(a) 2(y + e 2 ) 5 / 2 (y + e 2 ) 3 



k 2 
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(y + e 2 ) 3 



K 2 g + 



( e z + 1 - , 



(y + e 



2)5/2 (e 2 + !) 



H 



22 



(y + e 2 ) 5 ' 2 

j c 
(y + e 2 ) 2 



K + 



1 



(y + e 2 ) 5 ' 2 (y + e 2 ) 2 (e 2 + 1) 



£ 



We now express 0y in terms of 0fy by assuming the forcing 
to be incompressible, for simplicity. In that case, the following 
relations hold: /^(t) = R(r){f x + ef z ), and h 2 (j) - [<pT - e]f x + 
[y + e<pr]f z . Keeping only the terms at most proportional to e, 
we obtain 



(g + a 2 )Vn + 2e(g + a 2 )V, 3 + Oie 2 ) 



(B.8) 



Here, k - n/2 - T(a). To deri ve Eg. iBAi. we used f ormul a 
(2.173), (2.175) and (2.18) in Gradshtevn & Rvzhild dl965b . 
The results are written in terms of new integrals by making 
the substitution y = 2£(e 2 + 1)t 3 /3 and taking the limit f — > 0: 



+ a 2 )^) + 0{e 2 ) 



hi = (g + a )[benf/ u +^13] 

+e{-(a 2 + l)(g + a 2 )ifr n +1 
hi — b 2 a 2 ip\\ + 2gbaif/\3 + g 2 iff33 + e{-2ba(a 2 + l)i//u 

-2[a 2 b 2 + g(l + a 2 )]if/ l3 - 2gbafr 3 3] + 0(e 2 ) , 



y 3\28? + l-)) Jo y 2 / 3 



exp 



1/3 



at(e 2 + 1) 

L 35/2 |gg / 2 ^ 2 + 1 ) \ 

. J (e 2 + 1)1 3 J 
_ T(l/3) / 3 \!/3 3 2/3 e0 



.,2/3 



+ 0(^ /3 ) 



;) - 



= Qo 



\2^e 2 + \) l 



(e 2 + 1) 



+ <9(£ 2/3 ) 



(e 2 + 1) 



£ 



l_ r dy 

3 Jf^o y 



exp 



■lntf) 



3 5/3 e0 / 2<f(e 2 
(e 2 + 1)1 3 



+ <9(£ I/3 ln(£)) = £„ + <9(f /3 ln(£)) 



1/3 1 



We simplify iBAi by keeping terms up to at most second 
order in e. To this order, we have: 



y 



1 +b z 



R(a) = (y+ a 2 ) 
7r/2- arctan) — — ) 



2eba + e 2 a 2 



+ C^e 2 ) = k - -= 



7 + a A 



For simplicity, we assume that (pij(k y ) = (pij(-ky) and conse- 
quently neglect all the terms proportional to an odd power of 
b only (but keep those proportional to ba). Recalling that £ is 
a small parameter, we will keep only the terms proportional to 
£ -1 / 3 and £° and the leading order term in the e expansion. The 
dominant contributions then give us: 



V11 = - 



1 



V12 



V22 = - : 



2g(g + a 2 ) 
2b 2 „ 

^^5/2*0^° + ■ 

3 2 'W 



2g 



K 2 - 



2u4 



2*0 



2e 2 b 



(B.5)By using J B . 8b and Eq. dB.7i for the turbulent viscosities and 
by keeping only the dominant term for each component of the 
forcing, we obtain equations dl4> given in the main text. 
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0(e 2 ) 



(B.7) 



ba 



2eb 4 



H17 — 



2g 2 (g + a 2 ) 
2b 2 „ 

3 2 'hb 2 



g 5 ' 2 (g + a 2 ) 



KQ0Q + ■ ■ 



Here, g — 1 + b 2 = k 2 H jk 2 and the dots stand for higher order 
terms in £ and e compared to those which were kept. 



